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We show how to expand the free energy of a matrix model coupled to arbitrary matter 
in powers of the matter coupling constant. Concentrating on v uncoupled Ising models — 
which have central charge v/2 — we work out the expansion to sixth order for v = 1, 2, and 
3. Analyzing the series by the ratio method, we exhibit the spin-ordering phase transition. 
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1. Introduction 

One of the unsolved questions in matrix models Q is, how do you introduce matter? 
This is an important question. If the matrix model is viewed as two-dimensional euclidean 
quantum gravity, the matter is the statistical system which is coupled to the random 
surface. If the matrix model is viewed as bosonic string theory, the matter represents the 
target space in which the string propagates. And at least one phase of nonabelian gauge 
theory can be reduced to a matrix model with matter || . 

A general framework for coupling matter to matrix models has not been found. The 
only exactly solved cases are the open chain of matrices, which gives unitary models with 
central charge c < 1, and the one-dimensional case. The region c > 1 remains almost 
completely unknown, although it is easy to formulate matrix models with any central 
charge. Recently there have been a number of attempts to develop alternate schemes that 
one hopes could deal with matter in this region @,|J. The method and results of are 
promising, but one is nervous about modeling random surfaces out of about ten squares or 
triangles, which are rather unwieldy, especially if the surfaces are to be branched polymers. 

In the study of conventional spin models, low and high temperature expansions have 
proven valuable. Maybe such expansions can prove useful for spin systems coupled to 
lattices as well. Before Mehta solved matrix chains ||, Itzykson and Zuber [|] suggested 
expanding the two-matrix model in low temperature series. This suggestion has been 
forgotten, but the problem of coupling c > 1 matter to a fluctuating surface remains 
unsolved. Here we develop a low temperature expansion for matrix models. In each order 
n in the matter coupling constant, we get arbitrarily large surfaces in which all the links join 
equal spins, except for n "bridge" links. These bridges connect the several "blobs," open 
surfaces on which all spins are frozen equal; two blobs connected by one or more bridges 
have unequal spins. This is the random surface analog of the ordinary low temperature 
expansion. There are two advantages of this expansion: one can expand a system with 
arbitrary central charge to rather high (or infinite) order in the cosmological constant, 
obtaining surfaces which are — hopefully — close to the continuum; and one can use it to 
study the c — > oo limit. 
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2. The expansion 

The partition function of a matrix model coupled to arbitrary matter can be written 



(2.1) 



where d>j are hermitian N x N matrices. We consider cubic mo dels, with V(<j>) = <p 2 /2 + 
g<f) 3 / V~N- The type of matter is encoded by the number of matrices in the model and the 
matter coupling matrix Q, which depends on a = e - ^, the matter coupling constant. 

Here we specialize to v uncoupled Ising models. They are "uncoupled" only in the 
bare action, though; each one interacts with the fluctuating surface and they thus interact 
among themselves. At the critical point of such a model, the central charge of the matter 
is u/2. For v Ising models, we need 2 V matrices to represent every combination of spins at 
each site. The coupling matrix Q can be considered as the connection matrix of the target 
space graph; the graph is a //-dimensional hypercube, in which the connection strength 
between any two vertices is a raised to the power of the dimension of the lowest simplex 
which contains both vertices. For one or two Ising models, for example, the matrices can 
be written"!": 
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(2.2) 



In the spherical limit, the free energy of these models is 

F^(g,a) 



(2.3) 



(which we normalize for later convenience). The problem is to expand F^ v ' in powers of 
a. We begin by expanding = + z^a + z^a 2 + • • ■; of course, we will generate 
both disconnected and connected terms, but the former will be canceled by the logarithm. 
Consider z^: 

4 1) = (Tr0 1 2 ) = (0g)(0^), (2.4) 



where in the last expression the averages are with respect to a single matrix model. By 
expressing the traces in terms of components, we have reduced the calculation to two 
surfaces (blobs), each of which has uniform spins, with one link (bridge) joining the two. 



t A superscript in parentheses will refer to v. 
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Similarly, zi will contain contributions from two blobs joined by two bridges, or three blobs 
joined in an open chain. Obviously, we can repeat this procedure for all multi-matrix model 
averages, yielding contractions of one-matrix model Green's functions. 

The object of interest, then, is the one-matrix model average tensor (0^* •••</>^ ri ). 
Because of the (j) i— > U^<j)U symmetry, it can only depend on 8^, and we must keep 
separate the upper and lower indices In general, 

(«•••«)= E v^-x-w 

^en„ (2.5) 

where U n is the set of permutations of n objects. Because the matrix components commute, 
many of the coefficients in (^75|) are equal. We can characterize this by a mapping from 
permutations to the partitions of the integer n, / : 11^ — > P ni where f{n) is the set of the 
lengths of the cycles of permutation 7r; if /(tt) = f(n'), then = X n '. Therefore we define 

E «<*> [(^v.:SrW+2]£::SrW + -- : 

P eP n J (2.6) 

fM = /(tt 2 ) = ■■■ =p 

For n = 3, for example, we have 



Vi 0/? 2 ^3 / -^.ldft % % + «3,2 I V + d 03 dp 2 d Pi + dp 2 d pi dp 3 



(2.7) 



(for a given n, the k's will be numbered in lexicographic order of the partitions). 

The task now is to calculate the coefficients K n ,p{g)- (The reader may wonder why 
we do not simply use one-matrix model connected Green's functions |7j to represent the 
blobs. We wish that this were possible, but — despite one's intuition — blobs that make up 
a planar surface may themselves be counted by nonplanar Green's functions.) The most 
straightforward method is to contract ( |2.6|) with the various tensors T^"'p n (tt), obtaining 
a closed set of linear equations for K UyP , where the inhomogeneous terms are one-matrix 
averages of products of traces. This becomes quite cumbersome, though, when n gets large. 
Another method is to contract (|2.6|) with • • ■ A^™ , where A is some N x N hermitian 



tensor; this gives 

((TrA0) n )= /Un,^n, P TrA^TrA^--- {2.1 



pEPn 



where pi,p 2 , . . . are elements of the partition p, and \x n , v is the number of different per- 
mutations tc G fl n such that /(7r) = p. 

The averages on the left hand side of (|2.8| ) can be calculated by expanding the external 
field integral 

Z(g, A) = fv(f> e - T ^ v ^- A ^ (2.9) 



in powers of A. Fortunately, this integral has been computed in the spherical limit by 
Kazakov and Rostov || and by Gross and Newman |J using loop equations.^ Using Gross 
and Newman's notation, 

1 1 . 1 2 

F{g, A) = Jim — log ZG/, A) = - — 5 M^a + ^) " oil?-* ~ x ) + ^¥^^-3 



, 3 1 1 I O 

+ CTi (T_i + « / — ffj - - - log 3fir 



(2.10) 

where /x a = a/Aq, + x, cr fc = l/(A a + x) k / 2 , A a are the eigenvalues of A, and x satisfies 

the equation 

x = 1/(120) -y/fya^x). (2.11) 

We first expand x = xq + xiTr A + £21 (Tr A) 2 + a^Tr A 2 + • • ■; the correct root of 
fl2.11|) has xo = 1/(12(7) — 6g + • • •• Plugging x into ( 12.101) and comparing with fl2.8|) , we 
can calculate the required coefficients. There is only one subtlety: we do not want Green's 
functions where two external legs are connected directly to each other, as this would "short 
circuit" that surface; i.e., the number of dissident links would be one less than required. 
Since the Green's functions generated by T are connected, the problem only occurs in 
second order; therefore we must subtract the constant term from that Green's function 
before we exponentiate T to obtain Z. Proceeding confidently, we can now read off the 
coefficients. The first few are: 

^ = TekN + H 9 + ^ + N (212) 

1 1 



' This integral seems to be much harder for a quartic than for a cubic potential. That is the 
reason why we use a cubic potential in this work. 
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and xo satisfies xo = l/(12g) — y^3g/x. 

Equipped with the coefficients k, we can directly evaluate the free energies ( J2.3[ ) for 
any value of v. The a term is of course just 2 V times the free energy of the one-matrix 
model 0. In the Appendix, we give through order a 3 where the coefficients are 
given for arbitrary v as exact functions of v and g; in other words, we have included 
arbitrarily large surfaces. Beyond that, we have calculated F^ 2 \ and F^ through 

order a 6 , where the series coefficients have been expanded to order g 32 , meaning that we 
have included surfaces made of up to 32 triangles. We will present a brief analysis of the 
series in the next section, but first we will give a summary of the difficulties encountered 
in the calculation, and how we can check it. 

When one expands the interaction term in Z^ v > in powers of a, in each order one gets 
large numbers of graphs in target space, which is a zv-dimensional hypercube for v Ising 
models. The number of graphs increases rapidly with both v and the order in a. Many of 
the graphs are isomorphic, and so give the same results. The choice is between generating 
many labeled graphs, or many fewer unlabeled graphs. In the latter case (which is what 
we did up to order 3), however, one has to solve the difficult combinatorical problem of 
how many ways there are to embed each unlabeled graph in the hypercubic target space. 
The other major computational difficulty is in expanding the external field integral ( |2.9|) ; 
one could, perhaps, speed things up by dropping the quadratic term from V{<p) and then 



shifting 4> to induce it ||10|1. 

We can check FW, of course, against the exact result [TJJ]. For more than one Ising 
model, we can regroup the series into powers of g, and compare with small-surface expan- 
sions similar to the ones developed by Brezin and Hikami ||. This is a series to much 
higher order in a but much lower order in g than what we have. One may, of course, also 
use the low temperature expansion to check the small-surface expansion. 



3. Series Analysis 

To exhibit critical behavior, we perform a series analysis modeled after those in ||. 
Having expanded to sixth order, we do not expect to calculate the critical exponents, but 
we will learn something about the critical behavior. We first regroup the series into powers 
of g. The coefficient A n of g n is then the sum of the matter partition functions on all the 
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n-th order graphs, i.e., surfaces of n triangles; we know it only to order a , though. If we 
had calculated it to all orders in a, its asymptotic behavior at large n would be 



A n ^g- n n c (3.1) 



where the exponent £ = 7 s t r — 3. From the exact solution for v = 1 one finds |nj that 
((a) = —7/2, the pure gravity value, for all < a < 1 except the critical point a*^ = 
(2y/7 - l)/27 = 0.1589, which is where the spin ordering phase transition takes place and 
where £ = —10/3. This is how we will look for the phase transition. 

To calculate C( C X we wu l use a ratio method jl2H . We define TV = A n /A n -2, q-k = 
n(rn^ — rj?)/2, and for u > 1, 



\u— 11 / n i n\ \ u ~ 11 

nr v n 1 — (n — Zu + 2)r l n _ 2 1 



nq 



2p 

U ~ 1 ] On, _L 0\J U ~ 1 ] 



n 



(n-2u + 2)q [ :i 2 1] (3.2) 



2p- 2 

Ju] _ [u] I [u] 
°n ~ ' n / *in 

Naively, the asymptotic behavior as n — > oo should be 

rM « g~* [1 + {n- u )] 
S N « C [1 + 0(n- 



(3.3) 



This is true, however, only if there are no confluent singularities. Another difficulty with 
the method is that as u gets large, the coefficient of the n~ u term can get large as well. 

We start with one Ising model. In fig. [I] we plot sj^ ~ C(a) from the sixth-order 
series for various iterations of the ratio method, u = 3, . . . , 7. The known value of a*W 
is shown as a vertical line. At a = 0, £ is very close to its exact value for pure gravity, 
—7/2; for u = 4, for example, we have ( = —3.4991. These graphs are not exactly what 
we would have expected; the peak, for one, is much too high. Nonetheless, fig. |l| does give 
qualitative evidence for the spin-ordering phase transition. Moreover, as we increase it, 
we obtain increasingly accurate values for a* (the closest one being 0.1547 at u = 6),^ as 
well as lower peaks. The ratio approximants seem to deteriorate, however, after u = 4, 
developing a second peak. We will therefore use a compromise value, u = 4. 



t One can also get very good estimates for a* from the peak in the specific heat, C ~ 4-% log g* . 



In figs. 0-^, we plot £(")(a) for v = 1, 2, and 3, keeping terms up to orders 3-6 in a. 
The right-hand side of the graphs is given only for completeness, as there is no reason to 
trust the expansion when a is not small; if we had expanded to all orders in a, the curve 
would come back to —7/2 at a = 1, because at infinite temperature each spin fluctuates 
independently. Two things can be claimed with certainty: the plots are showing evidence 
for the spin-ordering phase transition which becomes stronger and more realistic as we 
expand to higher orders in a; and the critical point gets closer to zero as the central charge 
increases. 



4. The v — > oo limit 

As can be seen for the first three orders from the expression for the free energy in the 
Appendix, apart from an overall normalization, the coefficient of a n in the free energy is an 
n-th degree polynomial in v. This property, although not a priori obvious, can be shown 
to hold to all orders. The easiest way to see it is through the small-surface expansion @. 
There, the coefficient of g m in the free energy is a sum of the (one Ising model) partition 
functions of all surfaces of area m, with each partition function raised to the power v. Each 
of those partition functions is a 3m/2-degree polynomial in a, so its n-th power cannot 
have any higher power of v than u n . This proves the assertion. 

We can see immediately that the asymptotic dependence of T^ v ' on v becomes trivial 
as v (and therefore the central charge) approaches oo, since v will then simply be a multi- 
plicative renormalization for a. The critical exponents of the system will be independent 
of v, and the critical temperature will have the asymptotic behavior 

V 

This is already approximately true in our results for up to v = 3. In fig. [5] we plot o*^ 
for u = 1,2, and 3, showing the plausibility of the asymptotics (|4.1| ). It is also intriguing 
to note that in the third-order expression for the free energy in the Appendix, the graphs 
which appear as coefficients of v n in order n are all trees. If this behavior can be shown to 
hold in all orders, this will prove that surfaces in the c — > oo limit are branched polymers. 
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5. Discussion 



In this letter, we have shown how to expand a matrix model coupled to arbitrary 
matter in powers of the matter coupling constant. We have shown that already at order 
a 6 one can qualitatively observe the spin-ordering phase transition, with good quantitative 
results for the critical coupling. The immediate extension to this work would be to expand 
to higher orders (in a and g), and to refine the series analysis. The reader should know that 
the present work was carried out with rather simple algorithms programmed in the useful 
but very slow Mathematica language, running on an overloaded Iris 4D/480S; altogether, 
the calculations took about a day of real time to execute. For conventional spin models, 
the low temperature series are typically known to order 20 or 30, which gives accurate 
results for exponents; perhaps one could push the present expansion as far? 

Another possibility is to expand other thermodynamic quantities than the free energy. 
The magnetic susceptibility series, for instance, is known to converge rapidly. It is not hard 



to add a magnetic field to a matrix model One could, of course, also experiment with 
different types of matter such as Potts models. 

Finally, the limit v —>■ oo is worth studying. This limit becomes quite clear in the 
context of the low temperature expansion: one must calculate the v n coefficient of the a n 
term. Many of the calculations leading to the low temperature series drastically simplify 
in this limit. The combinatorics on a hypercube, for example, is much simpler when 
its dimension goes to infinity. There exist predictions for behavior of a random surface 
embedded in D — ► oo dimensions |T3[] which can be tested. In this way, perhaps we could 
learn whether the central charge c — > oo limit is universal. 



Appendix 

Here we give the free energy for arbitrary v to first three orders. ^ K n ,p,l is the coefficient 



As in eq. ( |2.3| ), the expression here does not include the homogeneous (pure gravity) term. 
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of the N £ term in K njP . 



1 



2" 2 '■' -~- 
1 

4 



1 r 
12 



(«?,!,! +2«; >1> iK 2i i_i + 2kJ )1) iK 2i 2 1 o) ^ 
(4,2,0 ~ 4,1 l) U ] tt2 + 

1,1, i (l + 6 (K2.1.-1 + K2,2,o) + 6 («2,1,-1 + «2,2,o)' 
+ («3,i,_| + 3k 3 , 2 ,_3 + 2« 3i3 _i )) z/ 3 + 

+ 3k 1, 1,1 (2«2,2,0«3,2,-f +4^2,2,0^3,3,-1 ~ 
- 2k x 1 iK2,l,-l - x i «2,2,o) ^ 2 + 

2 Kii+<3,-i)+ 3 +°( c4 ) 

The for n < were given in the text; the third-order coefficients are: 



«3,1 = 



= f 3g + 2x ^'dgx ^ N - 5/2 



V 8x 3 (3g - 4x 3 ) J 

+ 3k 11 ik 2 ,i,-iN- 1 / 2 + k* 1 N 3 / 2 

AT" 3 / 2 



«3,2 = 



«3,3 = 



16xq 
1 



+ «l,li«2,2, iV 1/2 



1 



32x5 8x 0y /3gx 



N -i/2 
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Figure Captions 

Fig. 1. £(a) for v = 1 (one Ising model), different iterations of the ratio method: u = 
3, . . . , 7. The vertical line marks the known critical point. 

Fig. 2. C(a) for v = 1, to orders a 3 , a 4 , a 5 , and a 6 

Fig. 3. Same as fig. z/ = 2 

Fig. 4. Same as fig. f = 3 

Fig. 5. The critical point a*^ for v = 1,2, and 3 
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